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Abstract 

The Thirring model and various generalizations of it are analyzed 
in detail. The four-Fermi interaction modifies the equation of state. 
Chemical potentials and twisted boundary conditions both result in 
complex fermionic determinants which are analyzed. The non-minimal 
coupling to gravity does deform the conformal algebra which in par- 
ticular contains the minimal models. We compute the central charges, 
conformal weights and finite size effects. 

For the gauged model we derive the partition functions and the ex- 
plicit expression for the chiral condensate at finite temperature and 
curvature. The Bosonization in compact curved space-times is also 
investigated. 



1 Introduction 



The response of physical systems to a change of external conditions is of eminent 
importance in physics. In particular the dependence of expectation values on tem- 
perature, the particle density, the space region, the imposed boundary conditions or 
external fields has been widely studied Nevertheless, many properties of such 
systems are poorly understood. The massless Thirring model 0, which is among the 
simplest interacting field theories, has already led to considerable confusion about its 
thermodynamic properties in the literature ^ ^. The reason is two- fold: Firstly, 
the computation of the fermionic determinant in the presence of a chemical potential 
and/or non-trivial boundary conditions is delicate, because the eigenvalues of the 
Dirac operator are generically complex. In section 3.1 we propose a regularization 
scheme via analytic continuation. We argue that the so-obtained determinant, which 
differs from previous results 0, leads to the correct equation of state. 
The second complication originates in the infrared-sector. An elegant infrared regu- 
larization, which is particularly well suited for the study of thermodynamic properties, 
is to quantize the model on a torus. Harmonic contributions to the current arise then 
naturally and taking them into account turns out to be crucial for a correct quanti- 
zation. In particular the so-obtained results differ from those gotten earlier using 
bosonization. This is explained in section 3.2. 

On another front there has been much effort to quantize self-interacting field 
theories in a background gravitational field For example, one is interested whether 
a black hole still emits thermal radiation when self-interaction is included. Due to 
general arguments by Gibbons and Perry this question is intimately connected with 
the universality of the second law of thermodynamics. The Thirring model (including 
the gauged version of it) is still solvable in curved space-time and we can study its 
properties in a background gravitational field. This provides us in particular with an 
elegant approach to the study of its conformal structure: Correlation functions with 
current- and stress-tensor insertions, which are gotten by functional differentiation 
with respect to the gauge- and gravitational fields, contain the necessary information 
to characterize the underlying symmetry algebras. To familiarize the reader with our 
approach we first rederive the conformal structure of the original Thirring model in 
section 3.3. We then show how a non-minimal coupling to gravity leads in a natural 
way to a modification of the conformal structure. In particular, very much as for a 
free scalar field the central charge in the fermionic formulation of the Thirring model 
is not unique. Furthermore, we find that the equivalence between finite size scaling 
and central charge of the Virasoro algebra holds only for a particular treatment of 
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the zero-mode sector in which a charge at infinity is generated automatically. This 
charge combines in a non-trivial way with the Weyl-anomaly of the determinant of 
the fiuctuation operators to reconcile the equivalence of the finite size scaling and the 
central charge. For certain values of the non-minimal coupling we obtain minimal 
models from interacting fermions. This is the subject of sections 3.4 and 3.5. 

The gauged Thirring model, which contains the Schwinger model {QED2) as a 
particular limit, is no longer conformally invariant but has a mass gap: The 'photon' 
acquires a mass = e^/(7r -|- ^g"^) via the Schwinger mechanism. It possesses a 
non-trivial vacuum structure which promotes it to an attractive toy model to mimic 
the complex vacuum structure in 4-d gauge theories. From our experience with the 
Schwinger model which is supposed to share certain aspects with one-fiavour QCD 
[y], we expect that gauge fields with winding numbers are responsible for the non- 
vanishing chiral condensate and in particular its temperature dependence. Configura- 
tions with windings, so called instantons, exist only for finite volumes and minimize 
the Euclidean action. They lead to chirality violating vacuum expectation values. 
For example, a non-zero chiral condensate develops which only for high temperature 
and large curvature vanishes exponentially. 

Since for particular choices of the coupling constants the model reduces to well 
known and well studied exactly soluble models there are many earlier works which are 
related to ours. Some of them concentrated more on the gauge sector and investigated 



the renormalization of the electric charge by the four-Fermi interaction |]T0| or the 



non-trivial vacuum structure in the Schwinger model P, |TT|. Others concentrated on 
the un-gauged conformal sector. Freedman and Pilch calculated the partition function 
of the un-gauged Thirring model on arbitrary Riemann surfaces [Q. We do not agree 
with their result and in particular show that there is no holomorphic factorization 
for general fermionic boundary conditions. Also we deviate from Destri and deVega 
[§] which investigated the un-gauged model on the cylinder with twisted boundary 
conditions. We comment on these discrepancies in section 3.2. 

Section 2 contains introductory material and in particular the classical structure of 
the model. 

Other papers are dealing with different aspects of certain limiting cases of the 
model considered here. In particular in H], the thermodynamics of the Thirring model 



has been studied and the Hawking radiation has been derived in |jT2[ . The equivalence 



of the massive Thirring model and the Sine-Gordon model in curved space has been 
shown in |T3|. Partition functions for scalar fields with twisted boundary conditions 



have been computed in [|14| and more recently in ||15 
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2 Classical Theory 



The gauged Thirring model in curved space-time has the Lagrangian 

where the gamma-matrices in curved space-time are related to the ones in Minkowski 
space-time as "i^ = e^al"' i V^ = c?^+ ioj^ — leA^ is the coordinate- and gauge covariant 
derivative and F^i, is the electromagnetic field strength. The gravitational field g^j,i, 
(or rather the 2-bein e"^, since the theory contains fermions) is treated as classical 
background field, whereas the 'photons' and fermions ip will be quantized. 
The classical theory is invariant under U{1) gauge- and axial transformations and 
correspondingly possesses conserved vector and axial-vector currents 



/ and j«'^ = ^7^73^ = ry'^,j-. (2) 

Here rj^^ = y/—g denotes the totally antisymmetric tensor. In fact, the conservation 
laws together with the relation @) between the vector- and axial currents imply that 
the currents are free fields 



2 -s/x 







(3) 



which is the reason that accounts for the solubility of the model 
presence of gauge- and gravitational fields 



even in the 



Of course, for any gauge invariant regu- 
larization the axial current possesses an anomalous divergence in the quantized model. 
Thus the normal Ua{^) Ward identities in the un-gauged Thirring model |]10| become 
anomalous when the fermions couple to a gauge field. 

The solution of the equations of motion is most easily presented by introducing 
auxiliary scalar- and pseudo-scalar fields, in terms of which the action takes the form 



S 



+g^-{d^ct>d,ct> + d^\dA) 



(4) 



Note that for later use we have allowed for different couplings of the fermionic currents 
to the scalar- and pseudo-scalar auxiliary fields A and 0, respectively. The original 
Thirring model is recovered for gi=g2=g, since then 
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= -\f^ + ^V^T^V^V^ + gfB^ + g^'B.B,, = d,\ - r//9,0 (5) 

is classically and quantum- mechanically equivalent to (|Tp, after elimination of the 
multiplier field B^. 

By decomposing the gauge field similarly as the S^-field as 

Af, = d^a - Vf.^d^if, so that Fqi = y^V^f, (6) 

and choosing isothermal coordinates for which g^j^i^ = e^'^rj^y, the generalized Dirac 
operator reads 

P = ^iP-h^G-la a ^-^F-^^,G+\. ^^^^^ 

(7) 

F = g^X + e a and G = g24> + e (p . 
Hence, if ^'0(2^) solves the free Dirac equation in flat Minkowski space time, then 

^(x) = e*^+*^^^-^>o (8) 

solves the Dirac equation of the interacting theory on curved space-time. The vector 
currents are related as 

The same relation holds for the axial vector current. 

Diffeomorphism invariance then leads covariantly conserved energy-moment um 
tensor 

Applying the variational identities in Appendix A one obtains after a lengthy but 
straightforward computation 

T^" = ^g^'^'F^'PF^p- F^^F^^" +^-[i)-i^^'D''^-{D^^'i))YH] 

+2V''(/)V"(/. - ^^"V>Va0 + (0^A) (10) 



where we have introduced the symmetrization A^^B"^ = ^(A^B" + A^B^). The first 
two hnes are just the energy momentum of the electromagnetic field, charged fermions 
and free neutral (pseudo-) scalars. The remaining terms reflect the interaction be- 
tween the fermionic and auxiliary fields. On shell T^^ is conserved as required by 
general covariance. Using the field equations for -0 and A its trace reads 

= -\f'^'F^, . (11) 
In particular for A^ = it vanishes, and the theory becomes Weyl-invariant. 



Symplectic structure: In the presence of both fermions and bosons it is convenient 



to exploit the graded Poisson structure [T^ 
{A{x),B{y)]^Y.jd^\ 



^f A{x) 5 5B{y) A{x) 6 5 B{y) 
50{z) S-Koiz) ^ 5t:o{z) 50{z) 



The sum is over all fundamental fields 0{x) in the theory . The sign is minus if one 
or both of the fields A and B are bosonic (even) and it is plus if both are fermionic 
(odd) fields. The momentum densities vro(a;) conjugate to the 0-fields are given by 
functional left-derivatives. A simple calculation yields the following momenta 

TT^ = 71^ = g^Jo + 2(9o0 and tta = g^jo + 2do\. 

In the following sections we are lead to consider the Euclidean version of the 
model. Then one must replace the Lorentzian 'y'^jg^u and by their Euclidean 
counterparts. For example, with our conventions the relation (Q) becomes 

r = -^V^f 

and the generalized Dirac operator in Euclidean space-time becomes 

p = e-^^e^^ ^ e^, where / = -iF + -f^G + (12) 

(see (|^) for the definition of F and G), instead of (^. Also, to recover the Euclidean 
Thirring model as particular limit of (§) we must set g^ = g1 = g'^- 
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3 Thermodynamic- and conformal properties 



In this section we analyze the quantum theory corresponding to the classical action 
@ without gauge fields, in fiat and curved space-time. The gauged model is then 
considered in the next section. Here we calculate the partition function, ground state 
energy, equation of state and determine the conformal structure of the un-gauged 
model. 

To allow for a non-vanishing U(l)-charge we couple this conserved charge to a 
chemical potential /z. For the finite temperature model the imaginary time must vary 
from zero to the inverse temperature j3 and the bosonic and fermionic fields must 
obey periodic and anti-periodic boundary conditions, respectively. We enclose the 
system in a spatial box with length L to avoid infrared divergences. 

We shall determine the dependence of the partition function and correlators on 
the metric. This provides us with an alternative approach to the conformal structure 
and its relation to finite size-effects. Also, it enables us to study the effect of non- 
minimal coupling to gravity in section 3.3. Hence we allow for an arbitrary metric or 
2-bein e^a with Euclidean signature. We can choose (quasi) isothermal coordinates 
and a Lorentz frame such that 

2a _ 2a f n\ ^ 2a 

QtJiu = e g^^ = e ^j, ^ = e tq, 

where r = ri + itq is the Teichmueller parameter and a the gravitational Liouville 
field. Space-time is then a square of length L and has volume V = Jq d'^x^/g. We 
allow for the general twisted boundary conditions for the fermions 



V 



ij{x',x' + L) = -e^-^^^+'^^^^Vla:",^^). ^^^^ 

The parameters and (3i represent vectorial and chiral twists, respectively. We could 
allow for twisted boundary conditions for the (pseudo) scalars as well fl^ e.g. 



(l){x^+nL,x^+mL) = (f){x^,x^) +27r(m+n). However, to recover the Thirring model 
for equal couplings we must assume that these fields are periodic. For cr = 0, T = i[3 / L 
in which case V = jSL, and for ao = jSo = the partition function has the usual 
thermodynamical interpretation. Its logarithm is proportional to the free energy at 
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temperature T = l/ (3. 



3.1 Fermionic Generating Functional 

Twisted boundary conditions as in require some care in the fermionic path in- 
tegral. The subtleties are not related to the unavoidable ultra-violet divergences but 
to the transition from Minkowski- to Euclidean space-time. To see that more clearly 
let denote the space of fermionic fields in Minkowski space-time with chirality ±1. 
Since both the commutation relations and the action do not connect and S~ we 
can consistently impose different boundary conditions on and S~ . On the other 
hand, in the Euclidean path-integral for the generating functional 

ZfIv^v] = j V^^V^ JV9^'m+SV9{-n^+^'v)^ ^^^^ 
the Dirac operator 

^ \D+ 

exchanges the two chiral components of ip, i.e. p : S^. Thus, in contrast to 

the situation in Minkowski space the two chiral sectors are related in the action. Of 
course, the eigenvalue problem for ip is then not well defined. This is the origin of 
the ambiguity in the definition of the determinant. It is related to the ambiguities 



one encounters when one quantizes chiral fermions |T8[. Here we reformulate this 
problem in such a way that the determinant with chiral twists [jS ^ 0) can be obtained 
by analytic continuation. The resulting determinants do not factorize into (anti-) 
holomorphic pieces. In appendix B we give further arguments in favour of our result 
by calculating the determinants in a different way. 

Let us now study the generating functional for fermions in an external gravita- 
tional and auxiliary field. For that we observe that on the torus we must add a 
harmonic piece to the auxiliary fields to which the fermionic current couples in (^. 
More precisely, in the Hodge-decomposition of in contains a harmonic piece, 

= d^X - + ^/i^ with V^h^ = hi^,,] = 0. (16) 

More generally, allowing for arbitrary couplings of the various terms in (|TBp to the 
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fermionic current, we are led to add a term 



to the action (p. Note, that in isothermal coordinates, for which the metric has the 
form (p!3D, the harmonics are constant. The constant /i^ couple to the harmonic 
part of the current and are needed to recover the Thirring model in the limit gl = g^ = 
gl- Also, we shall see that the harmonic degrees of freedom are essential to obtain 
the correct thermodynamic potential. 

Finally we introduce a chemical potential for the conserved U{1) charge. In the 
Euclidean functional approach this is equivalent to coupling the fermions to a constant 
imaginary gauge potential Aq [0. 

As a consequence of the above observations the scaling formula (|l^ (recall, that 
F = g^\ and G = g2(f) when the electromagnetic interaction is switched off) is modified 
to 



p = e ^"^e^^^e^, where f = -igA + l5g24> + ^cr 

P = + i(b^ - "^[goh^ + Ai^]) and /i^ = -^^At 5^0- 



(17) 



This scaling property will enable us to relate the fermionic determinants and Green's 
functions of p and p. The spin connection ci) in ([T7| ) vanishes for our choice of 
the reference zweibein. The dependence of p on the chemical potential fi and the 
constant harmonic field /i^ cannot be gotten by the anomaly equation pO[. It must be 
computed by direct methods. For this we expand the fermionic field in a orthonormal 
basis of the Hilbert space 



n n 



where an,bn,an, bn are independent Grassmann variables. A basis is given by 

i^n±{x) = e*(P"'^) e±, where (p^)i = ^(^ + ± A + rii), (19) 
and the e± are the eigenvectors of 75. Recall that and represent the vectorial- and 
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chiral twists (|1^) respectively. The ipn+ and ipn- must obey the and S~ boundary 
conditions, respectively. These boundary conditions fix the admissible momenta 
in (JT^I). Since the Dirac operator maps into the Xn± must then obey the same 
boundary conditions as the 4'n:f Thus Xn±{x) is obtained from ipn±{x) by exchanging 
and p~ . It follows then that 



with 



a; 



2n r 1 1 n \ 

— - r(- + ai+ (Ji + ni) - {- + ao + (Jo + no) 

r /I o N /I 

— T ^(o + oi - /i/i + ni) - (- + ao - Po + '^o) 

Tq-L l 2 z 



(20) 



(21) 



Here we have introduced = — Qohfj, — fi^. To continue we recast the infinite 
product for the determinant in the form 



1 



n,. 



(22) 



where g^'^ is the inverse of the reference metric 



and 



fl/x + ''Z^'^' with [f]^ = - 



1 

1"0 



1 



\t\ 

-n 



(23) 



The logarithm of the product (|2^) can in turn be written as the derivative at zero 
argument of a generalized zeta function. Indeed one easily verifies that for 



n 

we have (formally) 



(24) 



(25) 



However ({s) is divergent for s < 1. These divergences can be regularized as follows: 
We compute ({s) for s > 1 and subsequently define its value for s < 1 by analytic 
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continuation. 

Assume for the moment that is real or equivalently that there are no chiral twists 
f3fj_ and chemical potential /i. Then ({s) has a well defined analytic continuation to 
s<l via a Poisson resummation Indeed, writing ({s) as a Mellin transform 



(26) 



the generalized Poisson resummation formula 



(27) 



applied to the integrand in (|26|) yields after integration over t 



(2^ 



The zero mode with ra^ = is eliminated because for s > 1 it does not contribute. 
After this analytic continuation ({s) and C{s) are now regular at s = 0. More precisely 
C(0) = and 



C'(0) 



TT 



-1 



VaY^ {df^un'^n'') 2 exp[-27rm''cp 



\v{r)\ 



-Cl 

Co 



(0,r)e[ '^M(0,r) 

L Co J 



(29) 



Here we made use of det[c{ip)'^] = det[{ipy], which follows from ^(0)=0. 
For complex the Poisson resummation is not applicable and C'(0) cannot be cal- 
culated by direct means. To circumvent these difficulties we note that the infinite 
sum ( p4D defining the (^-function for s > 1 is a mereomorphic function in c. Thus we 
may first continue to s < 1 for real and then continue the result to complex values. 
Using the transformation properties of theta functions the resulting determinant can 
be written as 



-ai + Pi 
o-o — Po 



(0,r)e 



-ai - Pi 
ao + Po 



(0,r). (30) 



This is the main result of this section. 
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It can be shown that this determinant is gauge invariant, i.e. invariant under 

+ but not invariant under chiral transformations, — >• + as expected. 
Furthermore, it transforms covariantly under modular transformations r — r + 1 and 
r — s> — 1/r. In other words, det ip is invariant under modular transformations if at 
the same time the boundary conditions are transformed accordingly. The exponential 
prefactor is needed for modular covariance and is not present in the literature [Q. It 
correlates the two chiral sectors and will have important consequences. In Appendix 
B we confirm (^) with operator methods. 

The last step in the calculation of the fermionic generating functional is the inclu- 
sion of the local contributions to the auxiliary- and metric field, i.e. the dependence 
of the determinant on A, and a. For this we introduce the one-parameter family of 
Dirac operators 

= Ce^^^^e^^. (31) 

9r 



We take the r-dependence of the metric as gr = e'^'^'^g. With / as defined in (JT^, this 
family interpolates between p and p. The determinant of the full Dirac operator is 



then obtained by integrating the corresponding anomaly equation |2^ : 

Sl , g' 



detip = det{ip) exp + |^ | y^^A^) , (32) 



where 



Sl = / ^g\na - aAa] (33) 



is the Liouville action. In deriving this result we assumed that / ^/gX = 0. This 
constraint on the zero-mode of A (and similarly of 0) will be discussed below. Actually, 
for our reference metric the Ricci scalar TZ vanishes and the Liouville action simplifies 
to — / \^aAa. However, the above formulae hold for arbitrary reference metrics 
and arbitrary Riemannian surfaces. Furthermore, as expected for a gauge-invariant 
regularization, the function A and thus the longitudinal part of does not appear 
in the determinant. 

To complete the calculation of the generating functional we need to know the 
fermionic Green-functions S. Using the scaling property of the Dirac operator, eq. 
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(|3TD, it is easy to see that in an arbitrary background field S is related to S by 



S{x,y) = e-^(") S{x,y)e-^'^y\ 



Together with the relation (|32|) and the explicit form (p9| , |30D for detip this yields 
the fermionic generating functional 



ZF[ri,v] 



-Cl 

Co 



Cl 



(O,r)0 ^ (0,r) 



Co 
-.2 



g-/ ,(.)5(.,.).(.) . ^ 1^ I ^^^^ 



(34) 



By using the scaling properties of the Ricci-scalar and Laplacian (see appendix A) 
the exponent can be written in a manifest diffeomorphism-invariant way as 



96 

Here we used that on the torus TZ integrates to zero. On the sphere or higher genus 
surfaces the last formula is modified. 

The Integration over the auxiliary fields then leads to the full generating functional 
of the Thirring model. It contains all information about the thermodynamic- and 
conformal properties. This is the subject of the next two sections. 



3.2 Thermodynamics of the Thirring Model 

In this chapter we derive the grand canonical potential, equation of state and ground 
state energy for the Thirring model. For this we need to compute the partition 
function 

Z = J ShV<pV\ ZF[r] = f] = ^] e~^^, (35) 

where Zp is the fermionic generating functional (^4]) and Sb the bosonic action 

Sb = {2nf^r''h^K - J v^(A AA + 0A(/.) . (35) 

As it stands the partition function is still ill-defined unless we constrain the zero- 
modes artificially introduced in the Hodge decomposition of i?^ in ([l6|) . The choice 
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of the constraints is restricted by the symmetries of the system. In particular transla- 
tion invariance (or rotation invariance on the sphere) and covariance under modular 
transformations of the torus are symmetries which me may want to preserve by the 
zero-mode constraint. The constraint measure 

J dhodhiV(j)VXS{^)6{X) ■■■ = J dhodhiVgcpVsX 4>=^J^(j) (37) 

(and similarly for A) satisfies these requirements (The normalization by the volume 
in the definition of is needed such that the constraints and hence the partition 
function are both dimensionless) . For example, one finds the dimensionless partition 
function 



N^ = jvct> 5{4>) 



det'2(-A) 

for free bosons, where the prime indicates the omission of the zero-eigenvalue. 



(38) 



Integration over the hcirnionics: There is no restriction on the harmonic parts 
of the auxiliary fields and the Gaussian integral yields 



^ d^he 

— oo 



-Cl 
Co 



e 



-Cl 
Co 



e 



(A) 



47r^l + ^2/27r' 



where 



e 



^^■^ 'y ^in{n+u)A(n+u)+2ni{n+u)w 



(39) 



is the theta function with characteristics 



u 



1^ Vcti -h iri^(3y) and w = f ] ("o + iVo (^u - A^o) 



(40) 



and covariance 



T \ Tig;TQ 
-f'^'" 



9l 



27r + gl \ -47r - g^ 



-47r - gl 
9l 



(41) 
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Integration over A and 0: The integral over A, subject to the 5-constraint in (p?]), 
merely contributes one inverse square-root of the primed determinant of — 2A to the 
partition function and so does the integration over 0. In fact, to obtain the partition 
function of the Thirring model we divide Z by the corresponding partition functions 
A/q of the free bosons, eq. (^). Using (^) and (|3^) we obtain 



1 



AAo |r/(r)|2\27r + ^?2 



2vr + ^?2 ^ 



u 



w 



(A)e 



(1/24,1+932 )Si 



(42) 



where we have also used the scaling formula for the primed determinant of A [^O], |23 

(43) 



, det'(— aA) >/ n , r 1 / i 

det'(-A) = " • = log « • [4;^ 



with p being the number of zero modes of the operator. On the torus / ai = and we 
find 

det'(-aA) = - det'(-A), 
a 



which produces the extra factor ^/2^^+~g^■ In the Thirring model limit (72 = 9o and 
the square-root in (^) disappears. 



Zero-temperature limit: To investigate the thermodynamics of the model we 
assume space-time to be fiat and that T = i[3/L. Then 

1 Z 

is the grand canonical potential. First we analyze the low temperature limit of Vt. For 
/i = this yields the ground state energy. We observe that for T = if3/ L the covariance 
matrix A in (^) simplifies to 



zttA = — 



Id 



9f> 



Atv 27r -L -2 



9S 



9o -4vr - gf, 
-4vr - g^ gl 



(44) 



and has eigenvalues 
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(45) 



with corresponding eigenvectors 

vi = {-l,l) and V2 = {l,l)- (46) 
Also, the fj tensor (see ^3]) and /io (see |l^ in (|40| ) simphfy to 

° ^/^^ and a - -^ ^ a 

For — oo the saddle point approximation to the Gaussian sum (|39|) defining the 
theta-function becomes exact and therefore using that 

log|r?(r)p — ' p-^oo 



we find 



/3 ^ oo = -777 - + ^ 



TT 

mm 

2L nez2 



27r + 47r /zL 



^ {^2-m-;7^;^(/3i + ^)}' (47) 
2-n 

for the zero-temperature grand potential of the un-gauged model. Here the chemical 
potential and chiral twist enter only through the combination /?i+/iL/27r. Let us now 
discuss the potential in the various limiting cases. 

i) No chiral twist, /9i = 0, and vanishing chemical potential: Then VL{(3 — > 
00) coincides with the ground state energy. The minimum in (^7]) is attained for 
ni = n2 = [\ + ai\ and we find 

TT 27r 2tx 1 2 

Eo(L,«i, A = 0) = + Y2^rT^^"^ " '2 + "'^^ ■ (48) 

Only for anti-periodic boundary conditions, that is for ai = 0, does this Casimir energy 
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coincide with the corresponding result for free fermions. For gl > Att the Casimir force 
is always attractive whereas for gf^ < 4n it can be attractive or repulsive, depending on 
the value of ai. The result (^) is in agreement with the literature |Q. For example, it 
coincides with De Vega's and Destri's result if we make the identification ujj^j^ = 2TTai 
and 1//?J3J3 = 1 + gl/2TT in formula (42) of that paper. 



ii) Small twists and chemical potential: For small (3i and /i the minimum 
is assumed for = and the potential simplifies to 

r^,n \ vr 27r 27r 2 

^^^^^^ = -6L+T2^"^ (49) 

and does not depend on the chemical potential. For vanishing g^, the minimum of 
is attained for 

^1 = [- + "1 - Pi - 7;— J and n2 = [- + ai+ (Ji + —\, 
Z zn Z Zn 

where [x] denotes the biggest integer which is smaller or equal to x. This then leads 
to the following zero temperature potential 

n= - ^-^(A + ^)^ 

6L L 271 ' 

+ + (50) 



27r ^2 27T- 
For n = i3i = this reduces to the Casimir energy for free fermions with left-right 



symmetric twists and agrees with the results in |^ . 

Note, however, that for /5i 7^ we disagree with ^ . The difference is due to the second 
term on the right in (0). Let us give two arguments in favour of our result: The 
discrepancy arises from the prefactor appearing in the fermionic determinant (0). As 
discussed earlier this prefactor implies the breakdown of holomorphic factorization, 
a property which has been presupposed in 0. One can show that our results can 
be reproduced by starting with massive fermions and taking the limit m —>■ (see 
appendix B). 

The second argument goes as follows: Suppose that /5i = ai = 0. Then (|50|) simplifies 
to 
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n(fi ^ = -gj - ^(^) + ^(^ - 12 + ^1) . (51) 

For massless fermions the Fermi energy is just /i and at T = all electron states with 
energies less then /i and all positron states with energies less then —fi are filled. The 
other states are empty. Since dQ/dfi is the expectation value of the electric charge 
in the presence of fi we conclude that it must jump if /i crosses an eigenvalue of the 
first quantized Dirac Hamiltonian h. For vanishing twists the eigenvalues of h are 
just En = {n — ^)7i/L. From (^) one sees by inspection that the electric charge 

(Q) = — = 2[- + ^] = 2n for < fi< E^+i 
dfi Z Zn 

indeed jumps at these values of fi. Further observe, that in the thermodynamic limit 
L —!■ 00 the density 

L ^ ~2n + gl2^' 
reduces for = to the standard result for free electrons. 



Equation of state: We wish to derive the equation of state for finite T in the 
infinite volume limit L ^ 00. This may be achieved by interchanging the roles 
played by L and [3. More precisely, using that 



e 



u 



w 



(A) = A/det(zA-i) e^™-" 6 ^ (iA^^) 



u 



we find in analogy with the low temperature limit, that for L 
given by 



00 the pressure is 



vr 

— - min 

213 n^Z^ 



27r 



{ni + 712 + 2/5o}^ 



+ 



^^{n,-n, + 2a, + 2i^y 



271- 
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Here the minimum of the real part has to be taken. Again the minimization arises 
from the saddle point approximation to the theta function which becomes exact when 
L oo. For small twists the minimum is assumed for nj = and then 

TT 27r 27r 2 
op p zvr + 271 

becomes independent on the chiral twist jSo. As we have interchanged the roles of 
the temporal and spatial twists this is consistent with the earlier result that for small 
twists Q is independent of In particular, for ao = 0, we find the following equation 
of state 

This result is consistent with the renormalization of the electric charge which is con- 
jugate to the chemical potential. It shows that the thermodynamics of the Thirring 
model is not just that of free fermions as has been claimed in [^. Indeed, the zero 
point pressure is multiplied by a factor 2tt/ (27r + g"^). This modification arises from 
the coupling of the current to the harmonic fields. It is missed if only the local 
part of the auxiliary field is considered, which is the case if one quantizes the model 
in Minkowski space and then replaces the fco-integral in the Green functions by the 
Matsubara sum. This remark should also be taken seriously in four dimensions! Fur- 
thermore, we see that the 'pressure' p is real only for ao = 0, which is consistent with 
the finite temperature boundary conditions^]. 



3.3 Conformal structure 

In the first part of this section we derive the Kac-Moody and Virasoro algebras of 
the model (^) without gauge-interaction and prepare the ground for an extension, 
containing in particular the minimal models, in the second part. 
Recall ( [TT| ) that for = the theory reduces to a conformal field theory on fiat 
Minkowski space-time. To continue it is convenient to introduce adapted light cone 
coordinates = ± x^ and the chiral components of the Dirac spinor ip± = 
|(1 ± 75)'?/'. Then after substituting the classical equations of motion 

-'^This can also be observed in the Hamiltonian formalism IM. 
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2 (53) 
+id-{gi\ + 5f20)vr^+V^+ 

depends only on x~ and is therefore the chiral Noether current. Evaluating the 

Poisson bracket of the symmetry generator Tf = J dx~ f{x~ )T with the different 

fields yields the classical structure 

5f<p = fd.ct> ; 5/A = /9„A 
<5;^+ = /9_V^+ + ^^+9_/ ; 5fi^l = {5fi^^)' (54) 
5f3~ = fd^3~ +J-d-f ; 5/r__ = /9_T__ + 2r__9_/. 



Short Distance Expansions: Let us now determine the quantum corrections to 
these classical results. These are computed within the Euclidean functional approach 
from the short- distance expansions of the relevant n— point functions. We need not 
postulate Kac- Moody and Virasoro algebras in advance as has been done in [T^, pB[ . 
These structures are derive here. When comparing the classical with the quantum 
results one should keep in mind that the roles of ipl and ipl are interchanged when 
one switches from Minkowski to Euchdean space-time. In coordinates adapted to the 
holomorphic structure of the torus 

X = ifx° + ix^, so that = — {d^o — rdxi), 

the Dirac operator and the corresponding Greens function take the form 

where ^ = x — y, and the chiral components of the energy momentum tensor and 
current are given by 

T.. = ^(rT- + r-) = ^^T^'^ and j, = l(r/-/). 
From the conformal Ward identities 
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(0(xi) ■ • ■ 50{xi) ■ ■ ■ 0{xn)) = -f dz{Oix,) ■ ■ ■ 0{Xn)T,.) (55) 

we obtain the central charges and conformal weights directly from the correlation 
functions. However, because on the flat torus the expectation value of T^^ is constant, 
we need to compute at least the 3-point function to read off the conformal weights. 
As in the classical theory (see (^) the symmetric energy momentum tensor measures 
the change of the effective action F = log Z under arbitrary variations of the metric. 
On the torus there are two independent contributions. One being due to variations of 
the modular parameter r and its conjugate r which depend implicitly on the metric. 
The other is due to the variations of terms which depend explicitly on the metric. 
Since the chiral component T^^ is gotten by contracting T^^ with dg^y/df it follows 
that 



When doing metric variations it is always understood that we take the fiat space- 
time limit afterwards. The r variation is constant and may be discarded in the 
short distance expansion. Thus to analyze the algebraic structure we can work on 
any Riemann surface. This is not true for the finite size effects, which are global 
properties. This aspect will be analyzed in section 3.4. 

For example, taking three metric variations of the curvature dependent part of log Z 



with Z from (^2]) we find the following short distance expansions for the three point 



correlation function 



(27r)3 {u — vY{u — zY{v — z)' 



Substituting this result into the Ward identity (0) we obtain the central charge and 
the conformal weight of the energy momentum tensor 



c = 1 and hr^^ = 2 . (56) 

Note that the the central charge as well as the conformal weight are independent of 
the couplings g^ and g^. 

The conformal weights of the fundamental fields are obtained by computing the 
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fermionic two point function with stress tensor insertion 



Since Z ~ exp[F(7^^)], its metric variation vanishes after the flat space-time hmit 
has been taken. The variation of Sij can be found in appendix A. This yields 

1 , 1 _2 1 '^T^al 



h-ihn = h it = — I g. 



^1 2 I6n' 167r27r + c/2 



(57) 



Thus we have reproduced the classical results supplemented by additional and g2 
dependent quantum corrections. In the Thirring model limit 92 = 91=9, these terms 



add up to give the known anomalous dimension appearing in the Thirring model p6 



Furthermore, from (|57|) we may derive a condition on the couplings (71,5^2 if we insist 
on unitarity, i.e. on h > 0. We find 

In particular for > v^27r the conformal weights are positive for any real g2- 
Next we determine the Kac-Moody algebra of the U{1) currents. To derive the corre- 
lation functions with current insertions we couple the fermions to an external vector 
field, that is consider the 'gauged' model without Maxwell term. For example. 



e2 



The effective action with external vector field is then obtained by shifting the auxiliary 
fields in (|^) as 



5'20 fl'20 + ev? , gA giX + ea, (59) 

where = rj^duip + d^a and we have neglected the harmonic contribution to the 
external vector field, because it does not contribute to the short distance expansion. 
The resulting effective action does not depend on a due to gauge-invariance. To relate 
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the variation w.r.t. to that w.r.t. (p we use 

d,<P = V;Al, where A^ = A, - V A, 

is the transverse part of A^. We obtain the following short distance expansion 

,. 11 1 

{Jxjy) ' 



2-n2-n + gl (x - yY ' 
We read off the value k of the central extension in the f/(l)-Kac- Moody algebra 



k 



27r 



2TX + g, 



(60) 



The precise (^a-dependence of k (which can of course be rescaled to unity by an 
appropriate redefinition of the current) is related to a finite renormalization of the 
electric charge in the gauged Thirring-model which we will discuss in section 4. 
Finally, from 

/ ■ ■ T \ ^ ^ ^ 

[jxJy zz) ~ -47^2 2T,^gl[x-zY{y-zY 

we obtain hj = 1. 

To see how the left and right Kac Moody currents act on the fermionic fields we 
notice that after the integration over the auxiliary fields the A-dependence of the 
fermionic Green function factorizes as 

where g{x) = —ia{x)+'-y5Pip{x), jd = 2n/ {2TT+g'^) andm^ is the induced 'photon'-mass 
(see(P^)). Variation w.r.t. the A— field yields, after some algebraic manipulations, 
the U{1) charges 

<^^o = -^a + ^^) and ^-.0 = - ^^TT^)- (61) 

We have used the convention where the electric charge q+q is unity. In the Thirring 
model limit we can compare ( |6TD with the results obtained in For that we need 



to rescale the currents such that the central extension (|60D of the Kac-Moody algebra 



becomes unity jz Jl + gl/2Ti . It is then easy to see that we agree with Furlan 



et al. if we make the identification gp^ = g'^/4:TrJl + gl/2'K. 
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Non-Minimal Coupling: In section 3.1 we have analyzed the fermionic determi- 
nant in the presence of twisted boundary conditions. One may ask what happens if 
we introduce a local twist instead, that is 

ij{x)^'^{x) ; ^^(a;)"^e'"^(^), (62) 

which should be interpreted as a modification of the charge neutrality condition. The 
computation of the fermionic determinant in the presence of such twists is similar to 



that for a Weyl rescaling of the background metric (|3ll - |32D . Integrating the corre- 
sponding anomaly equation we find 

We will come back to the relation between the above determinant and charges at 
infinity at the end of this section. For the moment we use the analogy merely as a 
motivation to study the extension of the Thirring model obtained by coupling the 
A-field non-minimally to the background geometry. That is we consider the model 
@ again without gauge-interaction but with an extra coupling 



93 J n\. 



Then T in (^) is modified, 

r__ — > f__ = T__ + ^g,dl\. 

The corresponding modification of the classical conformal transformations ( |5^ ) gen- 
erated by the modified generator Tf = J dx~ f{x~)T are 

Sf(f) = 6f(f) , 6fX = 6fX - Y^-/ 

I - I (64) 

Whereas (p and ■?/'+ remain primary fields, A does not. This is in fact needed for 
consistency. Indeed, since if) is not a scalar under conformal transformations generated 
by Ty, the term ~ / ip'^pip in the action is only conformally invariant if A transforms 
inhomogeneously like a spin connection. 

It may be surprising that the new symmetry transformations depend on the cou- 
pling constant g^ which is not present in the fiat space time Lagrangian. However, 
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the same happens for example in 4 dimensions, if one couples a scalar field confor- 
mally, that is non-minimally, to gravity. Although the Lagrangian for the minimally 
and conformally coupled particles are the same on Minkowski space-time, their en- 
ergy momentum tensors are not. The same happens for the conformally invariant 
non abelian Toda theories which admit several energy momentum tensors and hence 



several conformal structures p7| 



The current still transforms as a primary with weight 1, but the energy momentum 
tensor acquires a classical central charge 

Sff^^ = /a_f__ + 2f__9_/ - g'.dlf. (65) 

The corresponding commutators in the quantized theory with non-minimal coupling 
to gravity are calculated as explained for the minimally-coupled model. One finds 
that the quantum corrections to ( p^ ) are identical to those of the minimally coupled 
model and thus are gs ^ 0-independent. 

To summarize, we have obtained the following Virasoro x Kac-Moody structure: 
Central charge: 

c = 1 + 2^gl^x and Ut^^ = 2 (66) 
Kac-Moody level and charges: 

. h-1 



27r + g. 



1/1 27r , _ 1,^ 2tt . 



2' 2-n^gl' 27t + g. 



Conformal weights: 



^° " 2"^16^^^~T6^27+^"^"^V 



~ 167r^^ lQ-K2Ti + gl 2 ~ ^ ' 

Here some comment about unitarity is in order. It can be shown that with respect 
to the standard scalar product reflect ion-positivity holds for any real g^ ||29|| . 
However with respect to this inner product the Virasoro generators are not selfadjoint. 
Choosing an alternative scalar product |14] for which they are selfadjoint, positivity 
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does not hold in general for g^,^ Q. We give a more detailed discussion about unitary 
subspaces in section 3.5. 



3.4 Finite size effects 

When quantizing a conformal field theory on a space-time with finite volume one in- 
troduces a length scale. The presence of this length scale in turn breaks the conformal 
invariance and gives rise to finite size effects. It has been conjectured that the 



finite size effects on a Riemann surface are proportional to the central charge. For 
example, when one stretches space time, — >■ ax", then the change of the effective 
action is proportional to c: 

Tax -r^ = -gloga-x, (68) 



where x is the Euler number of the Euclidean space time. In this conjecture has 



been proven for a wide class of conformal field theories on spaces with boundaries. 
The only important assumption has been that the regularization respects general 
covariance. In this subsection we shall see that the equivalence does hold only for a 
particular zero-mode treatment, which differs from (0). 

The only global conformal transformations on the torus are translations which do 
not give rise to finite size effects. Also, the Euler number vanishes and according 



to (pq) the finite size effects are insensitive to the value of c. For that reason we 
quantize the un-gauged model (H) on the sphere where the global conformal group is 
the Moebius group. 

An effective method to compute finite size effects has been developed in It 
is based on the following observation: Any conformal transformation z w{z) is a 
composition of a diffeomorphism (defined by the same w) and a compensating Weyl 
transformation g^^y — >• e^'^g^y with 

e- = ^^, z = ^ + ^x^. 
dz dz 

Therefore, choosing a diffeomorphism invariant regularization one has 

= ST^iff = 6Tconf — S^Weyl- 

The change of the effective action under Weyl rescaling is 

/P(A0) detiipg) exp{-SB[g]) 



ST Weyl = - log 



/D(A0) deti^p^exp{-SB 
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where Sb is the bosonic action (pG]). Since on the sphere there are no harmonic vector 
fields the term ~ h"^ in Sb is not present. Imposing the conditions (pT]) we obtain 



, y Sl gl f ,^ Svr, , det'A 
*r„.„ = log - - — + I / - -) + log (69) 

To evaluate ( |69D one introduces the 1-parametric family of Laplacians 

A, = e-2""A 

interpolating between A and A. Integrating the corresponding anomaly equation 
we end up with 



Consider now a dilatation w{z) =az. Then, the conformal angle is constant, (T = loga, 
and {IZ — Stt/I^) = 0. Then the first term in (^) vanishes and the finite size effect 
does not depend on gl. It is given by 

3 



"^r = -^iog« / \/^^ = - logo 



and does not agree with ( p5D since c in ( pBD depends on g^. On other Riemannian 
surfaces one would find the same result. Note that the finite size scaling comes from 
the middle term ~ log a J \fgR in (|70| ) which is topological in nature, while the short- 
distance behaviour of the energy-momentum correlators is controlled by the remaining 
two terms in (|70| ) which are insensitive to the topology. In that sense finite size scaling 



and the central charge are complementary. There is a way to match the two results 
by adding the term 



2 

-| / v^7^A7^ 



to the effective action. With this new effective action the short distance expansion 
of the energy-momentum correlators does not depend on g^ any more and the cor- 
responding central charge equals that obtained from the finite size scaling. However 
such a term would correspond to a non-local counter term to be added to the regu- 
larized action. 
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3.5 Charge neutrality and unitary subspaces 

In this subsection we show how the equivalence between the central charge and finite 
size scaling can be restored, provided the partition function is replaced by an average 
over un-normalized expectation values of charges at infinity. In fact it turns out that 
the J T^A-term, ie. the non-minimal coupling to gravity, itself can be given the 
interpretation of a charge at infinity if the zero-mode constraints (^) is replaced by 
a non-translation invariant sum over charges at infinity. 

The hint comes from inspecting the fermionic weights (|67D, which shows that 
iplx) and ipgs 

have the same conformal weights. We can therefore 
consistently put a charge at infinity with a corresponding modification of the charge 
neutrality condition. The non- vanishing two-point function is now {il!g.^{xyip{x)) . It's 
coincidence limit jg.^ is again a primary field with conformal weight hj = l. 
On the other hand, including a charge at infinity into the definition of the partition 
function we have 



3-5fl . gSTrgsAKo) 



Zg^ = fVs^VsX ZF[v = V = 0]e' ... 

^0 _ (71) 

= Zo exp[167r^^3G'o(^o,'^o)] (recall that V^cj) = ^{(fyVcl)). 



To continue we need to determine the coincidence limit of the scalar Greens function 
Go{x,y), i.e. to regularize the composite operator exp(aA) appearing in ([7I|). The 
normal ordering prescription 

oA(x) 



works well on the whole plane [^, On curved space we must be more careful 
when renormalizing this operator. The required wave function renormalization is not 
unique but it is very much restricted by the following requirements: First we take as 
reference system (the denominator in (|72D ) one with a minimal number of dynamical 
degrees of freedom since we do not want to loose information by our regularization. 
Second, the renormalized operator should have a well-defined infinite volume limit. 
Finally, the regularization should respect general covariance. These requirements 
then force us to take as reference system the infinite plane with metric gij,,^. The fiat 
metric 6f^u is not permitted since it leads to a ill-defined expression for (exp(aA)). 
With this choice the normal ordering in (ff^) is equivalent to replacing the massless 
Green function in (fflD by 
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Gl'^ix, y) := Golx, y) + ^ log [^l^s^{x, y)]. (73) 

Here s(x, y) denotes the geodesic distance between x and y. The occurrence of the 
arbitrary mass scale /i comes from the ambiguities in the required ultra-violet regu- 
larization. On the 2-sphere with constant Ricci scalar TZ we have 

Gl^%x,x) = -^[\og[^] + l]. 

The expectation value (: e*'^^^'^^^"-' :) then transforms under a constant rescaling z — > 

az as 

^. g8.93A(«o) ^ ^. g8M3A(?o) exp[87r(73' log(a)], (74) 
and therefore gives an extra contribution 

^r^s = ^ log(a)x, 

to the finite size scaling of the effective action. Adding this contribution to ( [fOD above 
we see that this is precisely the piece needed to restore equivalence with the central 
charge for any real or imaginary g^. 

More generally we can define the functional integral as an average over all possible 
charges at infinity: assume imaginary. The (un-normalized) expectation values are 
then given by 

nO„,(x.)) = - / VscpV\[-j= dk : e^'^^^?^) : ] l[OaX^^) e"^-- (75) 
i=i ' Z J V27r J J 

Here denotes the ?7(l)-charge of the operator Oj. In particular the partition 
function on is 



g-87rg3A(5o) . g-Ss[A'] 



where Aq is the zero mode and A' the excited modes of A(a;). The middle term in the 
above integrand is the zero-mode part of 5*^. The zero-mode integration yields a delta 
function (5(/c -|- iSirg-^i) and thus the g^ J 7?.A-term itself acquires the interpretation of 
a charge at infinity, due to the presence of the zero mode. The 'extra' charges e*'^'*''-^'^'' 
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assure the charge neutrahty of the partition function. For the general n-point function 
(|75[) the Aq- integration yields 



6{k + Sifig^ + ^a,), 

i=l 

where the sum of the ?7(l)-charges of the operators in ( [75|) enters. In particular, for 
neutral states, for which {J2 oii + iSvr^fg = 0), must be zero and no extra charge at 
infinity is present. 
Finally, using 

-^\dk e^^^(f») = 5(A(eo)), 
the averaging over all possible charges can also be written as 



PA5(A(eo)). (76) 

It is easy to verify that if the action has translation invariance in the target space, 
then the constraints ( [76l) and (|37D are equivalent and the correlation functions do not 
depend on the chosen base-point ^o- However, in the present case (|76D clearly breaks 
translation invariance (or rotation invariance on S"^) and the zero-modes constraints 
are inequivalent. Although we have assumed an imaginary g-^^ our results apply for 
any g,,. For particular values we recover the (unitary) minimal models, provided 
screening charges are included for the n-point function with n > 2. In particular 
for (yfg = I/a/IStt and gi=g2 = we obtain the Ising model with = = |. 



4 Gauged Thirring-like Models 

In this section we extend the model by gauging the global ?7(l)-symmetry. Contrary 
to what one might think, many aspects of the gauged model are actually simpler 
as compared to the ungauged model. In particular the thermodynamical properties 
are independent of external conditions like chemical potentials and twisted boundary 
conditions. The reason is that the model is closely related to the Schwinger model, for 
which the spectrum consist solely of a neutral, massive particle. On the other hand, 
the gauge interaction complicates the analysis, because the U{1)- bundle over the 
torus allows for gauge field configurations with winding number, so called instantons. 
These, in turn, imply fermionic zero-modes which trigger a chiral symmetry breaking 
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and therefore a non- vanishing condensate. This is the subject of the second part 
of this section. In the first part we discuss the partition function to which only 
topologically trivial configurations contribute. 

To see how the fermionic generating functional (|34D is modified, we decompose a 
general gauge potential on a torus as 



(77) 



where the last 3 terms correspond (as for the auxiliary field B^) to the Hodge de- 
composition of the single valued part of A in a given topological sector, that is the 
harmonic-, exact- and co-exact pieces. The role of the toron fields has recently been 
emphasized within the canonical approach . In the Hamiltonian formulation they 
are quantum mechanical degrees of freedom which are needed for an understanding 
of the infrared sector in gauge theories. Also, in ||36| it has been argued that the Z^- 



phases of hot pure Yang-Mills theories |^ should correspond to the same physical 
state if the toron fields are taken into account. The first term in (^) is an instanton 
potential which gives rise to a non-vanishing quantized flux. As noted above configu- 
rations with non- vanishing flux do not contribute to the partition function due to the 
associated fermionic zero modes. We can therefore assume Ai = for the moment. 
The fermionic generating functional is obtained from (^) by simply shifting 



Qohf, ^ etf, + goh^ = Hf, , g-,X^ea + gA = F and g2 



5f20 + e(p = G, 



which leads to 



27r(Vra'"'/3M/3<.-2i/3iao). 



1 



■ / r){x)S(x,y)ri{y) _ 



exp 



WW 
1 







-ai + I3i 



(O,r)0 



-ai - Pi 
ao + Po 



(0,r) 



(7J 



with a^ = a^ - H^ - /i^. 

To compute the partition function we must switch off the sources rj and f] in ( |75D 
so that 



Zo = J J dH(fhVipV(j)VX Zf[0,0] e 



(79) 



where now 
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Sb = {2Txf^gr''KK 

+ j V^(^v^AV - AAA - - g,R\ 



10) 



Note that we have kept the non-minimal couphng of the A-field to gravity as in 
section 3.3. Since Sb and the fermionic determinants are both gauge invariant and 
thus independent of the pure gauge mode a in ([77|), it is natural to change variables 
from Afj^ to This transformation is one to one, provided 



J Vd'f = J VdOi = and et^ G [0, 1]. 



^1) 



In contrast to the auxiliary harmonic fields /i^ in (|16|), the toron fields et^ and et^ + n^ 
with integer are to be identified, due to gauge invariance [§]. The measures are 
related as 



VAi, = JJ2 dtodtiVipVa, where J = (27r)2det'(-A). 



(82) 



In expectation values of gauge invariant and thus a-independent operators the a- 
integration cancels against the normalization. This simply expresses the fact that in 
QED the ghosts decouple in the Lorentz gauge. 

As we shall see shortly it is advantageous to integrate first over the toron fields. By 
using the series representation of the theta functions one computes 



d\et)Q\ '''l(0,r) e 
Co 



-ci 
Co 



(0,r) 



^3) 



Since the result appears always together with the T^-function factor in ( ^l)) it is con- 
venient to introduce 



/2^|r/(r) 



in the following expressions. The result (|83|) does not depend on the /i-field and hence 
the /i-integration in ( ffOD becomes Gaussian and yields a factor l/47r so that 
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Zo = 7rfi:det'(-A) e^^^^^'' J Vsiip(j)X) I V9GAG-Ss[h=o]^ ^^^^ 



where we inserted the exphcit expression (^2]) for the Jacobian. Now we see why we 
did well integrating over the toron fields first. It has washed out the dependence on 
the boundary conditions and chemical potential in (|83|). 

The integral over A, subject to the condition in (|37D , decouples completely apart from 
the non-minimal coupling to gravity which modifies the Liouville factor and yields 



one inverse square-root of the determinant of — 2A in (84). Thus 



Zo = fi:7ry'2Fdet'(-A) e(^3+i/24.)5^ 

■ J Vs{v<P) e^/v^^^^-^^[''=^=''l, ^^^^ 

where we have used (^31) . The ip -integration in contrast, leads to a finite renormaliza- 
tion of the dynamically generated 'photon' mass 



rriy J 

where 



n 271 + 

plays the same role as the r^'-mass in QCD The determinant obtained from the 
if — integration factorizes as 

det'(A2 - A) = det'(-A) ■ det'(-A + m^). 

This factorization property is not obvious since all determinants must be regulated. 
But it holds for commuting operators and in the zeta-function scheme. Then the 
partition function simplifies to 

Zo = ^^^^(det'(-A)det'(-A + m^))"^ exp {{gl + ^)S, 



We can go further by using the scaling formula for the determinant of A ||20[ and the 
known result for the determinant of A PT| which together yield 



det'^(-A) = roL|r^(r)|^^ exp ( - ^Sl). (87) 
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Thus we obtain the following partition function for the general model (^) on curved 
spaces: 



Zn = V2nV ; — 7—7- 1 exp ( ( h gf^) St 

Again we have factored out the partition function A/q for free auxiliary fields. The 
formula (|88|) shows explicitly that in the topologically trivial sector the theory is 
equivalent to a theory of free massless and massive bosons with mass m^, even in 
curved space-time |jl3[ . 



The appearance of in should be interpreted as renormalization of the electric 
charge induced by the interaction of the auxiliary fields with the fermions. After 
summing over all fermion-loops this leads to an effective coupling between the photons 
and the 0-field and in turn to a modified effective mass for the photons in (^). In the 
limit (72 this mass tends to the well-known Schwinger model result, — >■ ej ^Jv: 

We have already mentioned that the chemical potential coupled to the electric 
charge has completely disappeared from the partition function. This does not come 
as a surprise since the only particle in the gauged Thirring model is a neutral boson. 
This has no charge which may couple to the chemical potential. Also, if the partition 
function depended on \i then the expectation value of the charge would not vanish, in 
contradiction to the integrated Gauss law. The result (|88|) provides therefore another 
test for our result (|30D for the fermionic determinants with chemical potential. 
The final result is also independent of the chiral and non-chiral twists. The normal 
twists have been wiped out by the toron integration. In fact the chiral twists are 
equivalent to a chemical potential and therefore the above remarks concerning the 
chemical potential apply here as well. Did we assume holomorphic factorization for 
the fermionic determinant then the partition function would depend on the chiral 
twists. 

We conclude this subsection by giving the explicit formula for the partition function 
on the fiat torus. The zeta-function regularized massive determinant is expressed by 

det'(-A + m2)^ = ^e-5f'W, 

where 
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C'(o) = E 



7CL 



-.Ki{m-^L\J {n, n)) — 



(n, n 



4n 



^9) 



and (n, n) = Qijn^n^ is the inner product taken with the reference metric, and the 
sum is over all (n*) G Z2 with the origin excluded. For g^y = 5^^, in which case the 
partition function has the usual thermodynamical interpretation, the result reduces 
to one derived previously by Ambjorn and Wolfram 
infinity we recover a result in |19| 
simplifies then to 



In addition, if L approaches 
The free energy for ri = and on flat space 



f=-^logZ=^C'(0). 
with C'(0) from ( ^9]) and the particular choice for the parameters. 



4.1 Bosonization of the gauged Thirring model 

We pointed out in section 2 that for gi = g2 = g the classical theory (Q) reduces to the 
gauged Thirring model. The same is true for the quantized theory on the torus if in 
addition we set g^ = g. More precisely, the Hubbard- Stratonovich transform |^0| of 



the Thirring model is just the derivative coupling model (^) with identical couplings. 
In the process of showing that we shall arrive at the Bosonization formulae for the 
gauged Thirring model on the curved torus. We shall see that only the non-harmonic 
part of the fermion current can naively be bosonized and that for this part the rules 



of the un-gauged model on flat space time need just be covariantized. 

For that we calculate the partition function (|79D in a different order. First we 
integrate over the auxiliary fields. To understand the role of A and we introduce 
sources for them. Thus we study the generating functional for the correlators of the 
auxiliary fields 

Here 



is the action of the full theory, p is the Dirac operator in ( \[7[ j with all couplings 
set equal and Sb the bosonic action (pOj). Since A and integrate to zero (see (0)) 
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we may assume the same property to hold for the sources. The integration over the 
auxihary fields is Gaussian and yields 



9, A 



1 



where 



is the action of the gauged Thirring model on curved space-time and 

^0 



(90) 



(91) 



27rdet'('-A) 



(92) 



comes from the integration over the auxiliary fields. 

Let us first consider the partition function, that is set the sources to zero. Comparing 
( pOD with (|86|) and using (|87|) we easily find 



2 47r 



where 7 is the mean field (see (pT])) and we used ( p2D and 
neutral scalar field 7 is found to be 



e-i/^-^"''y'P7 5(7)e-^^ (93) 

. The action for the 



le 



1 



Since (|93|) holds for any (and thus for the non-harmonic part of any A^, because of 
gauge-invar iance) we read off the following bosonization rules: 



1 



i 1 



Jl + g^jl-K 



(94) 



where prime denotes the non-harmonic part of the currents. Thus, only the non- 
harmonic parts of the currents can be bosonized in terms of a single valued scalar 
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field. To bosonize their harmonic parts one would have to allow for a scalar field 7 
with winding numbers. On the infinite plane the harmonic part is not present and 
we may leave out the primes in (0). If we further assume space time to be fiat we 
recover the well-known bosonization rules in What we have shown then, is that 



for the gauged model on curved space time the bosonization rules are just the fiat 
ones properly covariantized and with the omission of the zero-modes. 

Since (|93D holds for any gauge field the current correlators in the Thirring model 
are correctly reproduced by the bosonization rules (Q). To see that more clearly we 
calculate the two-point functions of the auxiliary fields in the Thirring model ( pO] - 
p^ ). For that we differentiate ( pO] ) {ip is treated as external field) with respect to the 
sources and find 



1 

{(j){x)(j){y)) = -Go{x,y) + — I {Go{x,u)j^.^{u)Go{y,v)j^.^{v))T, 



(95) 



where Go is the free massless Green-function in curved space-time and the integrations 
are over the variables u and v with the invariant measure on the curved torus. Here 
(. . .)t are vacuum expectation values in the Thirring model ([91|). Alternatively we can 



calculate these expectation values from (84) and (|85|) , where the fermionic integration 
has been performed and find 

{\{x)\{y)) = ^Go{x,y) 

{mm) = ^Go{x,y) + ^{l-^)^{xMy). 
Comparing this with the result (^) we see at once that 

\Goix,u)rJu)Goiy,vmv))T = (96) 



2 

mz 



{Go{x,u)j^.^{u)Go{y,v)jl^{v))T = -^(m'^^ip{x)(p{y) - Go{x,y)y 

These correlators express the gauge invariance and the axial anomaly (j^.^) = —m^Aif 
in the gauged Thirring model. They can be correctly reproduced with the bosoniza- 
tion rules (|94D. They are not reproduced with the ones derived for the un-gauged 
model is! 



36 



4.2 Chiral condensate 



The chiral condensate is an order parameter for the chiral symmetry breaking. How- 
ever, on the torus its expectation value, whose temperature- and curvature depen- 
dence we will here compute would vanish if topologically non-trivial gauge field con- 
figurations were absent. There is a useful classification of the gauge configurations 
corresponding to the number of fermionic zero modes they give rise to. If we let 
A; = n+ — n_, where n± counts the number of zero-modes with positive/negative 
chirality, then we have 

/c = ^/ d?x^^ai{p\x) = ^d^xv^'F^,, = (97) 

which establishes a relation between the number of fermionic zero modes (or, more 
precisely the number of zero modes with positive chirality minus the number of those 
with negative chirality) and the first Chern character of the bundle. Also from ( P?]) 
one immediately concludes that the flux must be quantized in integer multiples of 27r. 
This is really a consequence of the single valuedness of the fermionic wave function 
(cocycle condition). 

Recalling the decomposition ( fTZl) of the gauge field we now concentrate on the 
first term which is the instanton potential giving rise to a non- vanishing quantized 
flux $. Since 2-dimensional gauge theories are not scale or Weyl invariant, as 4- 
dimensional ones are, the instantons on a conformally fiat space-time are not identical 



to the fiat ones [^, As representative in the fc-instanton sector we choose the, up 
to gauge transformations, unique absolute minimum of the Maxwell action in a given 
topological sector. It has field strength eE^ = ^^/V . The corresponding potential 
can be chosen as 

e4 = ei;.-$r//9,x, where ei^ = -^$(a;\ 0) (98) 

is the instanton potential on the fiat torus with the same flux but field strength 
y/g<^/V. The function x is then determined (up to a constant) by 

V9y - \j9j = v^Ax. (99) 
The solution of this equation is given by 
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Xix) = -li^e-'nix) = lj d'y^)Go{x,y)e~"^^y\ (loo) 

where 

Go{x,y) = {x\—-\y) = lOi 

is the Green-function for —A. In deriving ( |10CI| ) we have used that -^{^/V) =0 which 
follows from the spectral resolution ( |101| ) for the Green function in which the constant 
zero mode 0o = V W of A is missing. 

Our choice for the instanton potential (^) corresponds to a particular trivializa- 
tions of the f/(l)-bundle over the torus In other words, the gauge potentials and 
fermion fields at x^) and x^+L) are necessarily related by a nontrivial gauge 
transformation with winding numbers 

x^ + L)- A^(a;°, x^) = d^a{x) 

+L) = -e^""(^) g2«(ai+/3i75) ^(^3,0^ ^ly I J 

For the choice (^8D we find 

ea{x) = — — X . 
ij 

Note that A is still periodic in x° with period L and if) still obeys the first boundary 
condition in (|T^ . To calculate the fermionic zero modes we use the square of the 
Dirac operator 

= {"f* D%. ) = ^^^.v^S'"^" - \^ + h"'"--''' (103) 



In a pure instanton and harmonic background {(p = a = 0) on the flat torus (|103|) 
simplifies to 

= -g^^-b^b, - ^75. (104) 

In other words, is the same in the left- and right-handed sectors, up to the constant 
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/V . Furthermore this operator commutes with the time translations which leads 
to the following ansatz for the zero-modes 

where we have assumed > 0. The choice of Cp is dictated by the time-like boundary 
conditions in (|14D- Inserting this ansatz into the zero mode equation p'^Xp = yields 

where 2/ = a;^ + — -(cp — Hq). 

fv 

This is just the differential equation for the ground state of a generalized harmonic 
oscillator to which it reduces for t = itq. The solution is given by 



These functions do not obey the boundary condition ( |102| ), but the correct eigenmodes 
can be constructed as superpositions of them. For that we observe that 

so that the sums 

where p = 1, . . . , k, obey the boundary conditions and thus are the k required zero- 
modes. Indeed, since {ip)'^ in non-negative there are no zero modes with negative 
chirality because of ( p.04| ). With (|97|) we conclude then that there are exactly k zero 
modes with positive chirality. Modes with different p in ( |105| ) are orthogonal to each 
other and the overall factor normalized them to one, so that the system ( |105D forms 
an orthonormal basis of the zero-mode subspace. For A; < the zero-modes are the 
same if one replaces e+ by e_. 

To compute the fermionic determinant in a given topological sector we again introduce 
the one-parameter family of Dirac operators 
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Pr = ^T^e-^' Pe'^, p = r {d, + iCj^. - ^eij - ^ [^m + ^^^) , (106) 
which interpolates between p and p, similarly as in (^Tl). But now 

with F and G from (0), contains an instanton contribution. Also note that p contains 
the instanton part A^. To compute it's determinant we observe that the simple form 
(|104D of —p"^ allows one to reconstruct its spectrum completely [2^, |^: 

^2 _ f degeneracy = k 

" \2n^/V degeneracy = 2k. 

The corresponding determinant is |^ 

det't,^) = iff". (107) 

To relate the determinants of p to that of p we again integrate the anomaly 
equation, which now reads 

^^^^^ = /d^.v^(/(.)+/'(.) - - P,(,,m, (108) 

where, due to the fermionic zero-modes, the projector onto the zero-mode-subspace, 
Poix,pl) = , Krir) = (109) 

pr 

appeared. For the deformed operator p'^ the first Seeley-deWitt coefficient is 



1 



-R^ + 75rA-G + [(1 - T)^g ^ + -Jts- (110) 



Integrating w.r.t. r [^] one ends up with the following formula for the determinant 



in arbitrary background gravitational and gauge fields: 
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dei'lJp = det ^ dei'm exp (^ + ^[ ^GAG) 

(2k f $2 



;iii) 



In deriving this result we used that / y/gx = 0. 

Now we are ready to compute the chiral condensate {ip'^ P+ip) . Observing that 
the fermionic Green's function anti-commutes with 75 one sees at once that only 
configuration supporting one fermionic zero-mode with positive chirality contribute 
to the chiral condensate 

where 17+ = P+i]. Earlier we have seen that these are the gauge fields with fiux $ = 27r 
or instanton number k = l. Thus the condensate becomes 




(^tp^^) = -yJ^-^ I n-)4{^)Ux)^M-)e-'^^'='\ (112) 
where exp(. . .) stands for the exponentials in ( |111| ). First we integrate over the toron 



field t. The t-dependence enters only through the zero mode and more specifically ipQ 
in ( p.05| ) with p= 1. Using the series representation for the theta functions one finds 

dhijl{x)i)o{x) = ^. (113) 

Note that the result does not depend on the chemical potential similarly as in our 
calculation of the partition function. To continue we observe that the term / y/gG 
in exp(. . .) vanishes because of our conditions ( ^T]) and ( pTl) on the fields and 0. 
Furthermore Sgik = I] = Seik = + The remaining functional integrals are 
performed similarly as those leading to the partition function and we end up with the 
following formula for the condensate 

(^^tp^^) = /^|^(^)|2g-2.VeV+2./y/v/ix^e-2(.0+e^)(x)-<x{.)^ _ ^^^^^ 
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The expectation value is evaluated with 

1 2 2 

Seff = I v^[^^(A^ - - A)v9 - ^0A0 - ^0A^ 
A formal calculation of the resulting Gaussian integrals yields 

• exp exp - — ■ — -Go(x,x) , 



(115) 



where 

and Gjn; Gq are the massive and massless Green-functions. Here we encounter ultra- 
violet divergences since Go{x,y) is logarithmically divergent when x tends to y. To 
extract a finite answer we need to renormalize the operator exp(a0) as explained in 
section 3.5. This wave function renormalization is equivalent to the renormalization 
of the fermion field in the Thirring model and thus is very much expected [0, ^ 
The fiat Green's function on the torus 



1 , 1 



1 , 1 „,....o/^.. ,Te+e 



|2 



where ^ = x — y, possesses the logarithmic short distance singularity 

Go(x, y) = -- log - - log (4vrVo|r^(r)n + 0(0- (117) 

Furthermore 

Go{x, y) ~ Go(a;, y) + 2x{x) - i / y^x + 0(0- 

With the prescription explained in section 3.5 we find that on the fiat torus G]^^ has 
now the finite coincidence limit 
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log 



47rVo|r/(r) 



:il8) 



To determine the chiral condensate we also need to determine K{x, y) on the diagonal. 
In a first step we shall obtain it for the flat torus. Its a-dependence is then determined 
in a second step. For a = and r = itq the Green function K has been computed in 
[H]. The generalization to arbitrary r is found to be 



m'iK{x, x) 



1 



2myLTQ 



coth(?2?)+ ' 



+ ^{~\og\vi—)\' + FiL,T)-HiL, 

/TT V T 



(119) 



27r ^ ■ • r 

where we introduced the dimensionless constant a = Lm^\r\/2n and the functions 



F{L,t) 
H{L,t) 



n>0 



E 



„>o VnFTa 
We used the abbreviations 

1 



+ 



g— 27riz-|-(ri) ]^ g27ri2_(n) ^ 



(120) 



(121) 



Substituting (|119|) and ( p.l8| ) into ( |115| ) with cr = we obtain the following exact 
formula for the chiral condensate on the torus with fiat metric (jfj^iy-. 



1 /m^L\ 



L\ 



32 

) exp ( 



exp 



2tt ^ ^c^Ltq 

{FiL,T)-HiL,r 



^ coth ^ 



2 r 



(122) 



where we used that on the fiat torus x = and V = V. Furthermore, we identified /z 
with the natural mass scale of the theory. 

To extract the finite temperature behaviour of the chiral condensate we take r = 
if3/L where (3 = 1/T is the inverse temperature. In the thermodynamic limit L — > oo. 
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Then coth(. . .) ^ 1, H ^ and the expression for the chiral condensate simphfies 
to 



{i-'P.i-h = exp [ - + ^F]. (123) 



^2 



Using 



F{P) ^ 7 + log ^ + ^ for a ^ oo, 



where 7 = 0.57721 ... is the Euler number, we obtain the zero temperature limit 

(V'tp^V) = -52^^^/(2-+^i)exp(^-^7) for T ^ 0. (124) 

For temperatures large compared to the induced photon mass F vanishes. Thus we 
obtain the high temperature behaviour 

2 

(^tp^^)T = -T(^)^exp(-^T) for T ^ 00 (125) 

It is instructive to discuss the various hmiting cases. For all gi = 0, i.e. the Schwinger 
model limit, the exact result ( p.23| ) simplifies to 

{^^P^.p)r = -T e-^^^^(^) ^ { Zf,%n., tZIo, (126) 

where now = e^/vr is the induced photon mass in the Schwinger model. This 
formula for the temperature dependence of the chiral condensate in QED2 agrees 
with the earlier results in 

Next we wish to investigate how the self interaction of the fermions affect the 
breaking. For large coupling g2 and fixed temperature the exponent in ( |123| ) vanishes 
so that 

{iI)^P+iI))t ~ , for T fixed, g2 00. 



Hence, for very large current-current coupling the chiral condensate vanishes. Or 
in other words, the electromagnetic interaction which is responsible for the chiral 
condensate, is shielded by the pseudo scalar-fermion interaction. 
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For intermediate temperature and coupling g2 we must retreat to numerical evalua- 
tions of the sums defining the chiral condensate in (|123|) . The numerical results are 
depicted in Fig. 1 

How does the gravitational field affect the chiral condensate? To answer this 
question we need to know the massive Green's function, entering in (|116|) , for non- 
trivial gravitational fields (for simplicity we assume T = 0). Let us first consider a 
space with constant negative curvature. Then Gm-y has been computed explicitly in 



4J]. Here we only need its short distance expansion, given by 



- — {27 + log (^— ) + ^(- + «) + ^i- 



a 



Ois')}, (127) 



where a'^ = ^ + and tjj{z) is the Digamma function. Substituting (|127|) into ( |116| ) 
we end up with the exact formula for the chiral condensate for constant curvature 



7^=o ■ exp 



r TT 2- '-^^ -1 . ,4 



128) 



The asymptotic expansions for large- and small curvature are easily worked out in- 
serting the corresponding expansions for the Digamma function We find 



{tP^P+^P)n = iij^P+ij) 



11=0 ■ exp 



TT 



12e2 



:7^ 



for < 1 



and 



(129) 



{iP^P+^P)n = (^^P+^)7^=o 



•2m2- 



exp 



r vr 
4^ 



n 



7rm: 



4e2 



■7 



I7^ 



for —> 1.(130) 



Hence the chiral condensate decays exponentially for large curvature analogous to the 
high temperature behaviour. However, the pseudo-scalars do not suppress the effect 
of the curvature in contrast to ( |125| ). Comparing the exponentials in ( |130| ) to ( |125| ) 
we may define the curvature induced effective temperature as 

Teff = (131) 

In passing we note that if we compare the prefactors, rather than the exponentials, 
we would write 

-7^)5 



47rV2 ' 



(132) 
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The latter identification actually coincides (up to factor of 2) with the Hawking 
temperature of free scalars in de Sitter space 0. The correct identification involves 
the (dynamical) mass of the gauge field and is therefore not universal. From this 
observation we learn that the temperature associated with curvature depends on the 
matter content. Note finally that the non- minimal coupling (g^) has no effect on 
the chiral condensate. In Fig. 2 we have plotted the chiral condensate for arbitrary 
constant values of the curvature. 

For gravitational backgrounds with non-constant curvature we have to refer to 
perturbative methods for the calculation of the massive Green's function. Again 
we only need the short distance expansion of Gm^. For geodesic distances s small 
compared to the massive Green's function may be approximated by the Seeley- 
DeWitt expansion EH 



-| oo o 

Gmix,y) - ^^j:a,ix,y){ - 4'\ms), (133) 

where Hq is the Hankel function of the second kind and order zero. In particular 

i/f (z)-.-[log^+7] for z^O. 

171 Z 



Inserting (|133|) into (|116|) we end up with the following expansion for the chiral 



condensate in an arbitrary background 

{^^P^^U = (^tp^^)7e=o -exp - -(^)^5:a,(x)^^^-^ , (134) 

where we have used that a^ix) = 1. The first order contribution involves ai(x) = — 
and reproduces the asymptotic behaviour ([1291). Higher order contributions must be 



taken into account to uncover the effect of variable curvature. For this one has to 
substitute is the corresponding Seeley DeWitt coefficients into ( |134|) . These have 
been computed up to j = 5 . 



5 Conclusions 

In this paper we have elaborated on various features of the Thirring model as well as 
some of its extensions. In particular we found the dependence of the partition function 
on the chemical potential and the non-trivial boundary conditions for the fermions 
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on the torus. For that a careful analysis of fermionic determinants has been crucial. 
We have found that the familiar chiral anomaly of the UV-regularized two point 
function is also seen in the IR-sector as a breakdown of holomorphic factorization. 
This fact, which has not been properly taken into account previously, together with 
the presence of harmonic contributions to the current, leads to a modification of the 
equation of state due to the current-current interaction. We believe that our results 
could also be obtained in the bosonized theory, provided the usual bosonization rules 
are modified to include scalar fields with winding numbers, i.e. scalar fields with 
vahics in a compactified target space. 

Furthermore, we have deformed the conformal structure by allowing for different cou- 
plings in the transversal- and the longitudinal parts of the current-current interaction. 
This does not change the Virasoro- and Kac-Moody algebra, but modifies the con- 
formal weights of the primaries and in particular of the fermionic fields. Not all 
values of the coupling constants belong to physical theories, since positivity of the 
scalar product imposes certain restrictions on them. Our approach allows also for 
a non-minimal coupling of the longitudinal sector to gravity. While such a coupling 
may seem to be ad-hoc we gave some arguments that it might arise naturally when 
quantizing fermions in presence of a a background charge. We find that the central 
charge of the Virasoro algebra is sensitive to the non-minimal coupling. In particular 
c < 1 occurs for certain values of the coupling constant. However, we have not been 
able to derive constraints on this extra coupling without referring to the result by 
Friedan, Qiu and Shenker. We believe that an independent derivation of their result 
within a fermionic model would be most welcome. We have also established that the 
central charge controls the finite size effects only for a particular treatment of the 
zero-modes of the auxiliary fields which is equivalent to an average over charges at 
infinity. 

Finally we have considered the gauged Thirring model in curved space-time. We find 
that the partition function is independent of vectorial as well as chiral twists and the 
chemical potential. This result, which technically is due to the harmonic contributions 
to the gauge-fields, is in fact expected as a consequence of Gauss's law. Furthermore, 
using the (probably not so obvious) factorization property of the zeta-function reg- 
ularized determinants of commuting operators we find that the partition function 
can be expressed completely in terms of a single massive scalar field. The gauged 
Thirring model shows a chiral symmetry breaking which originates in the existence of 
fermionic zero-modes and thus in configurations with winding number (instantons). 
We have obtained explicit expressions for these instantons as well as the expectation 
value of the chiral condensate as a function of temperature and curvature. The con- 
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densate is exponentially suppressed for high temperatures and/or big curvature which 
is interpreted as an almost restoration of the chiral symmetry under these extreme 
conditions. Although temperature and curvature have qualitatively the same effect 
they cannot be identified. In particular the identification with the Hawking temper- 
ature for free scalar fields in de Sitter space does not hold in the present situation. It 
follows from general arguments that the chiral symmetry can not be restored for any 
finite temperature or curvature so an exponential suppression is most we can expect. 
In fact, it has been argued earlier, that the axial [/(l)-symmetry in 4 dimensional 
QCD also shows an almost restoration as a function of the temperature |^9|. Our 
results on the curvature dependence could motivate a corresponding investigation in 
QCD. Finally we note that the chiral condensate is linearly suppressed for large 
current-current couplings. 
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A Conventions and Variational Formulae 

Our conventions for the metric and curvature agree with those of Birrell and Davies 
[|]. We use the chiral representation 7^^^- = (Ti.Ym = for flat space with Lorentzian 
signature and 7^ = c^i^Tb = ""^2 in Euclidean space. Furthermore 75 = 75 = a^. 
In what follows we derive some variational formulae used in the text. Here denotes 
the space-time and Lorentz covariant derivative. 

Using the definition of the Christoffel symbols it is straightforward to show that 



1 

5^^^ = _Ye''a5e: ; = \{t''59,c. - g^'5g^p) (135) 

1 

2' 



= -g'-'^iDJgp, + D,6gp, - Dp6g,,) 



For some formulae related to the variation of the tetrad let us refer to ^ 
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1 (136) 
where t\ = - (e'^^^e.b - e\5e,^) . 

In addition we have 

5ujf,ab = -D;^^a6 - Ot^ab ] Oi^^ab = ]^(^''a(^^b{DJgpn " DpSQa,,). (137) 

When performing the variation of curvature dependent expressions we have used the 
identities 

g^'^6n,, = , where cj^ = g'^'^SF-^ - g^'^SF^^, 
and /^a;M„ = /V^{<7"'3v^A'^-V"A/^}55«^. ^ ^ 

Depending on the topology of space-time, the reference curvature TZ may be different 
from zero. In this case it is not possible to express the conformal angle a in terms of 
the curvature scalar. Nevertheless, to perform variations of cr-dependent expressions, 
the identity 

5{^n) = -2S{^Aa) (139) 

proves to be useful. 

Taking the variations of the equations 

^nG{x,y)^-5{x-y) and ^ ipS{x,y) ^ S{x - y) (140) 

for the scalar and fermionic Greens functions we may derive (up to contact terms) 
the following variational formulae 

^^^/ i-lg^'g'"' + g'''g^ndaG{x,u)dpG{u,y)^Sg,. 
SS^y (2S{x,u)YD^S{u,y)-D4S{x,uh%^r]'"'S{u,y)])^6g,.. 

Here all arguments and derivatives which are not made explicit in the integral refer 
to the coordinate u over which is integrated. Finally, we need the following formula 
for the variation of the inverse Laplacian 
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* i^f) = Z - *<^'Z^) - W I ^^"'^S^^f' (141) 

where V is the volume of space-time and / an arbitrary function. To prove this 
identity we note that for / G (KernA)-*- we have 



A 



Varying this equation yields 



A{6^f) = 6f-{6A)^f 



which may be inverted to give 



Varying the identity 



allows to replace the last term of (|142|) to obtain the required result (|141|). 



B Canonical Approach to the Partition Function 

In this appendix we compute the partition function for massive Dirac fermions in the 
canonical formalism. In the limit m ^ we confirm the result ([30| ) for the fermionic 
determinant with chemical potential in chapter 3. For massive fermions one cannot 
consistently impose chirally twisted boundary conditions. However, from the explicit 
eigenvalues (0) one sees at once that the chiral twist Pi and the chemical potential 
are equivalent. One can easily verify that this equivalence holds also for massless 
fermions in the canonical approach and that jSi ~ fiL/2TT. Let us therefore compute 
the partition function 

Z{(3) = Tr[e-^--^^-''^^--] (143) 

for massive Dirac fermions with chemical potential /i on a cylinder with (non chiral) 
twisted boundary conditions 
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^{x + L,t) = -e^2'™iV^(x, t). (144) 
For massive particles it is more convenient to use the Dirac representation 

7° = (^3 7^ = -i(^2, 7^ = 7°7^ = -(^1- (145) 

The Dirac field is expanded in terms of the eigenmodes of the first quantized Hamil- 
tonian 



idr —m 



(146) 



as 



^(a;,i) =I]^n,+&n + 5I^",-4, (147) 

n n 

where the iiJn,+ and ipn- are the positive and negative energy modes, 

_ ^-iuJnt-iXnX I _ iuJnt~iX„X 

Yn,+ "^71) Yn,— /I'^n; 

c„ = (2^„(^„ + m)L)-^ ("""aI ■ ^^^^^ 

The momenta A„ and frequencies Un are determined by the boundary condition (|144|) 
to be 

A„ = ^(n - ^ - ai) and Un = \Jni? + A^. (149) 

After normal ordering the 'positron' operators with respect to the Fock vacuum de- 
fined by H we find 

{H - /iQ) = Y^i^n - Mbr, + Y^i^n + /^)4c^n " Y.i^n + /^), (I50) 
n n n 

where the last c-number term represents the infinite vacuum contribution which must 
be regularized. To do that we employ the zeta function regularization. That is we 
define the zeta-function for s > 1 by the sum 

C(.) = EK + /i)-^ 
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which in turn defines an analytic function on the whole complex s-plane up to a simple 
pole at s = l. The analytic continuation is constructed by a Poisson resummation 

+ = — ^F(n), (151) 

n n 

where 



F(0 = e^-^«(^-i) J dy e'^y[^m^ + y^ + (152) 
and 771 = Lm, ji = Lix. Taking the Mellin transform of ( |152| ) we find 

F{C) = e2"^«(5-°i)_i_ J dy e'^y J dt f-^e~'^rh^W-t'^ 

= -^e2-^«(^"^) J dt f-'e-'^ ±Ko{fi^e + t') (153) 

2^ 27rjf(i-ai) f jj. J.S -til 



F diverges at ^ = since the Kelvin function Ki{z) ~ 1/z for small z. It follows that 
the n = term in (|151|) diverges. This divergence is regularized by subtracting the 
ground state energy of the infinite volume system. Indeed, because of the exponential 
decay of the Bessel function for large arguments, only the n = term contributes for 
infinite volume. So we find for the regularized sum 

EK + /.) =Yr^i:jdte he . (154) 



Now we perform the limit m — > 0. Only the most singular term in the expansion of 
the Bessel function contributes, hence 



(n2 + t2) 



^ E e2™(^--)v^n-i5_,_.^. (/in), (155) 



TT ' » 2 ' 2 



where Sa;b{z) is the Lommel function |5^. In particular for s = — 1 this function is 
5* = 1/z so that finally 
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TT 2n In 
^-^(«i-[«i + 2]) • 



Inserting this into ( |150| ) then yields the regularized expression 
V V 6L L ^ 



2^ 



For small n the normal ordering is /i-independent so that 

{0\:H-f^Q: |0) = + " [«i + ^])' = (0| : : |0) 

is independent of /i and coincides with the Casimir energy 

Let us now compute the partition function. Using ( p.58|) we easily find 



tr [e 

oo 

n 

n>[i+ai] 
oo 



/3:(i/-MQ): 



n (1 + g("-i+"i)e^'^) ■ 

n>-[i+ai] 



n (1 + g("-5+°i)e-'3^) 



n 



n>Ni-ai 



ra>- 



k(r) 



—ai 



(0,r) 



— «i 



(0,r) 



(156) 



(157) 



(158) 



(159) 



where we have used the product representation of the theta functions in the last 
identity and that q = e^'^*'^ = e~^'^^/^. A non- vanishing chiral twist Pi can now be 
included by shifting the chemical potential. Thus we have confirmed the formula ( PDD 
in the text. 

Note that for n ^ the zero-temperature limit of the grand potential is not equal 
to the vacuum expectation value of : if — fiQ : . For /i 7^ all states up to the fi- 
dependent Fermi energy are filled. For example, for ui < fi < 002 iii the limit (3 ^ 00, 
fl reduces to the expectation value of : H — fiQ : in the one-electron state. 
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Fig. 2: The chiral condensate as a function of the curvature and the 

couphng constant §2- 
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